Geometric multipartite entanglement measures 
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Within the framework of constructions for quantifying entanglement, we build a natural scenario 
for the assembly of multipartite entanglement measures based on Hopf bundle-like mappings ob- 
tained through Clifford algebra representations. Then, given the non-factorizability of an arbitrary 
two-qubit density matrix, we give an alternate quantity that allows the construction of two types of 
entanglement measures based on their arithmetical and geometrical averages over all pairs of qubits 
in a register of size A'^, and thus fully characterize its degree and type of entanglement. We find that 
such an arithmetical average is both additive and strongly super additive. 
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Quantum entanglement is arguably the most intrigu- 
ing feature of the quantum world, the hallmark of non- 
local correlations that draws the line between classical 
and quantum behaviour it is also the very resource 
that allows quantum computation 0] , and quantum pro- 
tocols such as teleportation and key distribution [1| 
to be performed. Since information can, in principle, be 
processed by quantum technologies via the manipulation 
of a given set of physical resources [l| , an important prob- 
lem to be addressed is the quantification and controlled 
manipulation of the degree of entanglement of a contin- 
uously interacting physical system [5|. In this context, 
we should be capable of performing a reliable measure 
in order to account for the quantification and processing 
of the system's degree of quantum correlations at a given 
time. Many efforts have been devoted towards this objec- 
tive, ranging from polynomial invariants density 
matrix properties @, I^, , and 
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In this Letter we give a geometric formulation of multi- 
partite entanglement measures. In so doing, we construct 
a Hopf bundle like mapping which allows us to obtain 
information about the entanglement of a pure quantum 
state. The idea of using Hopf fibrations has been explored 
by several authors in the context of the generalization of 
Bloch spheres to higher dimensions [l^ [3, [l7] . The two 
qubit case has an interesting interpretation in terms of 
the second Hopf fibration 17|. The three qubit case, al- 
though explored [15i, has not been fully understood. 
Here we give the formulation for the three qubit case and 
generalize it to an 7V-qubit system. This implies a tech- 
nical challenge since Hopf fibrations cannot be faithfully 
constructed in higher dimensions. We propose a general- 
ization in terms of Clifford algebra representations, and 
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build a Hopf-like map which gives information about the 
system's entanglement. Moreover, we introduce the per- 
mutation group to decode the information available to all 
qubits in the system. This geometrical picture provides a 
natural setting for the Meyer- Wallach entanglement mea- 
sure 0. We also propose, to avoid some inconveniences 
posed by the Meyer- Wallach measure, an alternate strat- 
egy to quantifying multipartite entanglement. We intro- 
duce the idea of probe quantities to characterize the non- 
factorizability of each two-qudit density matrix within a 
pure quantum state and build on it to construct two en- 
tanglement measures, that can also account for the case 
of mixed states, by means of their geometrical and arith- 
metical averages. We study their basic properties and 
implications. We find that the proposed arithmetical av- 
erage is both additive and strongly super additive. 



Two-qubit entanglement and the second Hopf 
fibration. — We construct the two-qubit case follow- 
ing Ref. but fixing the state norm to 1. The 
state vector can be written in the computational basis 
as 1^-) = ao |00) ai |01) -I- /3o |10) + /3i where 
ao,ai,/3o,/?i G C, and |aoP + \ai\^ + |/3oP + = 1- 
Although there has been a geometrical formulation of 
few-qubit entanglement in terms of Hopf fibrations in 
Refs. [l^ [l^, [l3|, this is limited to only three Hopf 
fibrations which would impose a limit on the system's 
scalability. Here we use a Hopf bundle-like mapping 
similar to the one reported in Ref. 15] but instead of 
using the geometrical picture of the spheres of the total 
space of the bundle we resort to the Clifford algebra 
representation. This construction allows us to avoid 
the limit imposed by Hopf fibrations and generalizes 
to higher dimensions. The relevance to a geometrical 
interpretation of states shall be dealt with later, as 
now we will emphasize its algebraic properties and 
build a map from the Clifford Cl{?>) representation 
(CT(3) - H ® H) to S"'. We introduce the pair of 
quaternions qi = ao -I- aii2, and q2 — + I3ii2- By 
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choosing the map tt : Cl{3) — > S' 



4. 



I*) = (<Zl,92) ^ (2(7291, kil' - k2|') = 

(2[/3oao + Piai] + 2[/3iao - /3oai]«2, 



|ao| 



l/3op-|/3i 



(1) 



we access ah the information available from the first qubit 
and a component related to the so-called concurrence C, 
a measure of the qubit system entanglement 5|| . Here x 
denotes complex conjugation. Thus, the map reads 



(2pi(oi) +Cii2, Pi(oo) - Pi(ii)) 1 



(2) 



where Pn(ij) denotes the one-qubit reduced density ma- 
trix of the n-th qubit. We point out that although we 
are keeping the map for its algebraic features, it is man- 
ifestly invariant under local SU(2) transformations and 
avoids the problem of defining charts as in the Hopf bun- 
dle picture. Such a map contains information about the 
first qubit entanglement, which should be symmetrical 
for the other qubit, and the issue here is how to recover 
it. In the two qubit case the other reduced density matrix 
can be recovered from the first one simply by exchanging 
ai with /3o- Of course this is equivalent to labeling the 
qubits in the other possible way, and in this sense it is 
a natural consequence of the arbitrariness introduced by 
the label we imposed. For higher dimensions, however, 
the number of possible labels is N\, so in principle we 
would have N\ different mappings. We now introduce 
the permutation group as a way of recovering the infor- 
mation for the other qubits. In fact, by permuting the 
particles we obtain the desired exchange, or equivalently, 
the other map: 

(2(aiao + /3i/3o) + 2(/3iao - /3oai)j2, |ao|^ + |/3o|^ - 

+ = (2P2(01) +C2j2, P2(00) - P2(ll)) ■ (3) 

As is to be expected, the term related to the concurrence 
is an invariant: the entanglement of the two particles re- 
mains the same, regardless of which one is measured. We 
will not go into the details of discussing the interpreta- 
tions of C (see e.g. Ref. [17i]), instead we shall focus on the 
fact that the map allows the construction of a measure 
of entanglement and that all possible maps are obtained 
through the permutation group. 

Three-qubit entanglement and the third Hopf 
fihration. — Although the three-qubit scenario has 
been discussed from the point of view of Hopf fibrations 
[1, 

[T^ . and also from the perspective of the twistor 
geometry formalism [YJ\ , there is not a proper systematic 
characterization regarding multipartite entanglement 
measures. Here we provide such a description. We 
use the permutation group to decode the information 
available from the Clifford algebra representation. The 
Hilbert space for the three-qubit system is the tensor 
product of the 1-qubit Hilbert spaces Tii ® 0.2 ® Ti-s 
with the direct product basis where a pure three-qubit 
state reads 1^-) = ao |000) + ai |001) + Po |010) -I- 



/3i|011) + JollOO) + 5i|101) + 7o|110) + 7i|lll), 
ao,ai,/3o,/3i,'5o,'5i,7o,7ie C, and jaop -I- |q;iP -I- |^oP + 
+ \So\^ + \Si\^ + |7oP + |7i|' - 1. We perform a 
similar parametrization to the third Hopf fibration but 
from the CZ(4) representation to S^: 

|*> = (01,02) = (gi -l-g2«4,g3 + 94*4); gi=Q;o + ai«2, 

92 = /3o + /3l«2, 93 <5o + <5l«2, 94=70+71*2, (4) 

which results in the map tt' : (5/(4) — > S^: 

^'(01,02) - (2020i,|0i|2-|02n= (5) 

(2Ci + 2C2i2 + 2C3H + 2C4Z6, |oip - |02n 



w_ithCi = ao<5o+<5iai+7o/^o+/3i7i, 
(/3i7o -JiPo), C3 = -PqSq + joao 



C2 — — ai(5o+(5iQ!o — 
- ("171 - (Si/3i), and 
(5i/3o + q;i7o- (/3i<5o-7iao), where |oip- |o2p = 
Pi(oo) -Pi(ii), and 1 = = |7r'(oi,02)p. Here Ci 

is clearly the off diagonal element of the reduced density 
matrix pi , but it is not trivial to see what exactly is the 
information that is available from the other C"s. In the 
two-qubit case, we had only one component, which was 
interpreted as the concurrence. In the three-qubit sce- 
nario there is certainly more information to look at. Con- 
sider the following qubit configurations: i) j^I'i) ® |^'23), 
C2 = 0, C3 = 0, Ci = 0, ii) l^-ia)® 1*3), C2 = 0, C3 ^ 0, 
Ci + 0, iu) 1*2) ® l^-w), C2 ^ 0, C3 = 0, C4 + 0, and 
iv) l^-i) ® 1*2) ® 1*3), C2 = C3 = C4 = 0. To some 
degree it may be tempting to say that for such configu- 
rations, the (123) case as we shall call it from now on, C2 
has information about the entanglement of (13) and C3 
of (12). An explicit calculation reads: i) case 1 (23): 
2Ci = 2pi(oi), 2C2 = 0, 2C3 = 0, 2C4 = 0, ii) case 
2® (13): 2Ci - 2pi(oi), 2C2 = 2a2C(i3) +262(7(13) , 2C3 = 
0, 2C4 = 4Im(a6C(i3)), "i) case 3 ® (12): 2Ci = 2pi(oi), 
2C2 = 0, 2C3 = 2a2C(i2) +262(7(12), 2C4 = 4Re(a5C(i2)), 
where C(i3) denotes the concurrence of the entangled two- 
qubit state (13). It is clear that C2 contains some infor- 
mation about the entanglement of (13) and Cj, about 
(12). We show next that, in contrast to the entangle- 
ment measure suggested in Ref. '^15], this information is 
not complete and that it is actually the = 3 permu- 
tation group which indicates the complete information 
available in the CI representation. We write down the 
explicit transformations induced by the elements of the 
permutation group of three elements as follows: 



a) (123) 

b) (123) 

c) (123) 

d) (123) 

(5o ^ /3o 
/3i -> 70 



(213) 
^ (321) 

(132) 
^ (312) 

70 /3i 
(5o ^ ai 



Pi 



5i, 5o ■ 

■ <5o. Pi ■ 
Po, Pq 

■ <5o, Po 



Po, Po ■ 
70, ^0 - 

■ ai, 61 

■ ai, Pi 
(231): ai - 



> So, Si - 
ai, 70 - 
' 70, 70 
* Si, Si 
■ Po, Po 



Pi, 
Pi, 
Si, 
70, 
So, 



e) (123) 

Si ^ Pi, Jo ^ Si. Hence, permuting 
the qubits is equivalent to redefining five new mappings 
through the equivalences obtained above: 



(123) -> (213): 

C'l = aoPo +_Piai + jqSq + Siji, C'2 = PiytQ - aiPo + 
jiSq - Sijo, C3 = joao - PoSo - iPiSi - 7iq;i), 
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-PoSi - 7oai + f3iSo + 71 ao, |oiP - I02I' 



P2(00) - P2(ll) 

(123) -> (321): 

C( = ao_ai+f3if3o+SiSo+^oji,_C2 ^ Siao-Soai+^i(3o- 
JoPi, C'3 = /3iao-/3oQ!i+5i7o-^o7ii (^4 = -70"! + 

(^O/?! - Pl^O + 71 ao, |oip - |02p = P3(00) - P3(ll) 

(123) ^ (132): 

C( = aoSo + 7o/3o + ^i^i + ^171, = 7oao - PoSq - 
if3iSi - 7iai), C3 = (5iQ!o - 6oai - (/3o7i ~ 7o/3i)i 
C4 = /Soi^i - 7oai - {(3iSo - 71 "o), |oip - |o2p = 
Pi(oo) - Pl(ll) 

(123) ^ (312): 

C{ =_Q;g/3o+/3iai+7o'5o+^i7i: C'i = 7oao-^o/3o+7iai- 
i5i/3i, C3_=_/3iao-ai/3o+7o<5i-i5o7ii = -/3o<5i- 
7oai +/?i(5o + 7iao, |oip - |o2p = P2(oo) - P2(ii) 

(123) ^ (231): 

C[ = agai+/3i/3o+<5i^o+7o7i7 C2 = /3i_Q;o-/3oai+7i^o- 
701^1, C'3_= (5iQ;o-(5oai+/3i7o-/3o7i, (^4 = /^o'^i - 
7oai - /3i(5o + 7iao, |oiP - |o2p = /93(oo) - P3(ii) • 

From this we see that (123) and (132) have information 
about the first qubit, (213) and (312) about the second 
qubit, and (321) and (231) about the third qubit. Al- 
though the density matrix information is the same for 
(123) and (132), the information of the entanglement 
available in (123) is not the same as that of (132) [2^ . 
and this is why the use of the complete permutation group 
of N — 3 particles is what gives the complete informa- 
tion available from the tripartite system. For higher di- 
mensions, and once a map such as Eq. ([5]) is constructed, 
permutations will always yield similar maps with the first 
and last components related to the reduced density ma- 
trices of the system. A natural way of extracting the in- 
formation hidden in the parametrization is given through 
the norm of the elements of the base space. We thus de- 
fine the quantity Ki as 



I2C3P + I2C4P 



Ki = I2C2I -I- |zo3| 

1- |2Ci|2- ||oi|2- |02| 

|2 |„ „ |2 



= 1 - 4|pi(oi)P - |pi(oo) - Pi(ii)l'' , 
which, after some manipulation, is equivalent to 
K, = 2(1 - Tr[p?]) . 



(6) 



(7) 



It is encouraging that this coincides with the Meyer- 
Wallach-Brennen quantities defined to measure entangle- 
ment in Ref. [23| , which can be seen by rewriting the mea- 
sure built by Meyer and Wallach through an analysis 
of invariant polynomials and the linear entropy. Indeed, 



an arithmetical average over all permutations yields the 
Meyer- Wallach-Brennen measure 



1 ^' 

^-l^E^^-l^E2(l"Tr[p^])- 



iV! 

1 1 

N 



N N 



(8) 



where pi are the reduced density matrices of one qubit 
accessible through the action of the permutation group. 
In the last equality we have used the fact that each Ki 
(fi) appears with multiplicity N — 1. For higher dimen- 
sions we can build similar maps to Eq. ^ using higher 
Clifford algebra representations in such a way that 



I2C2P + I2C3I 
l-|2pi(oi)P- 



+ |2C4p + ...- 
Ipi(oo) " Pi(ii)l 



I2C..-1H = 



(9) 



Multipartite qubit entanglement. — The case of dimension 
N > 3 requires special attention. It is known that the 
Meyer- Wallach-Brennen measure is not at all successful 
at accurately measuring the entanglement when iV > 3. 
The first qualitatively different scenario we find is the 

= 4 case jl^l, which allows the additional possibility 
of being factorized as a (2-|-2)-qubit state. In particular, 
the case of the direct product of two EPR states gives 
A/ = 1 , which appears as an unexpected result since this 
is a semi factorizable state. This is so because, for such 
states, the corresponding reduced density matrices are 
equal to ij and hence the information about the entan- 
glement available for each qubit is the same: maximal 
entanglement. 

Scott [1^ has addressed this issue by considering not 
only one qubit reduced density matrices, but the reduced 
density matrices oim < N qubits, corresponding to other 
bipartitions: 2 -I- 2, 3 -I- 2, etc. He has proposed the 
measure 



Qm(*) 



where CjY — 



E 



1 - Tr 



(10) 



\s\-- 



N 



The choice A^ = VN/2\ would ac- 



count for all possible bipartitions, and as any s-partition 
is necessarily included in a bipartition, then it should be 
enough just to consider bipartitions. In a similar spirit, 
Love et al. 0] proposed a geometrical average of the 
term over which the sum is performed in Eq. pop as a 
way of characterizing global entanglement. Nevertheless, 
these type of measures have the inconvenience that they 
do not yield one for generalized GHZ states, as all their 
reduced density matrices have purity equal to 1 /2 yield- 
ing a value equal to one only for the case of two qubit 
reduced density matrices. 

Although this can, in principle, be overcome by intro- 
ducing a proper normalization factor, we seek an alter- 
nate measure that can be averaged in a more satisfactory 
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way, namely, that yields one for GHZ states and gives 
a non-zero value for maximally mixed density matrices. 
Our proposal is that the arithmetical and geometrical 
averages of such a quantity are enough to characterize 
the entanglement of a quantum state, quantifying it and 
deciding whether the entanglement is global (among all 
parties) or not (bifactorizable, A''-factorizable). 

Based on the factorizability of a given density ma- 
trix, we propose an alternate strategy which is both eco- 
nomic and effective at quantifying multipartite entan- 
glement. The basic idea is to characterize the degree 
of entanglement of a system through a probe quantity, 
say 7-", that measures the degree of non- factorizability 
{pab = Pa® Pb) of each pair of qudits. 

Note that factorizability is only equivalent to separa- 
bility in the case of pure density matrices. For the case 
of mixed density matrices, a factorizable density matrix 
p = Pa®Pb can indeed be rewritten as a separable matrix 
X^KPa ® Pb^ however, the converse is not true (see e.g. 
the Werner states). We focus on qubit systems {d — 2). 

First, we have to find a suitable candidate for this task. 
Several candidates can be thought of, e.g., the mutual 
information [l^, the trace distance, and others 0]. For 
the sake of completeness in our analysis we consider two 
main quantities: i) The first quantity, which we term 
as the quasi- concurrence Qc{pab), is based in the same 
eigenvalues \i of the concurrence [14], but considers a 
different combination to that defined in Ref. 14]. Thus, 
we consider, in decreasing order, the eigenvalues of the 
matrix ^JpABPAB, where pab = {02 ® a-i)PAB{pi ® 0-2), 
in order to define Qc{pab) — Xi + X2 — X3 — X4. Note 
that Qc is non-negative, ranging from zero to one, and 
that for pure states is equivalent to the concurrence 
C{pab) = max{0, Ai — A2 — A3 — A4}, as we only have one 
non-zero eigenvalue. Also, Qc{pab) = for factorizable 
density matrices, since all their eigenvalues are the same; 
and all the reduced density matrices within a general- 
ized GHZ state \GHZ) = (1/^2) (JO)®^ + \lf^), yield 
Qc{pab) = 1, thus being a well suited quantity for our 
purposes, ii) The second quantity, Tr{A, B), has a more 
direct physical significance. This is proportional to the 
von Neumann's mutual information 

J^r{A,B) = ^(^S{pa) + S{pb)^S{pab)), (H) 

where S{p) is the von Neumann's entropy associated to 
the density matrix p Q. For example, for a generalized 
state \GHZ), all the two qubit reduced density matrices 
yield Tr{A,B) = 1/2; a W-state \W) = (l/v^)(1100) + 
JOIO) -I- JOOI)) has J'r{A,B) ~ 0.46 for all {A,B), and a 
fully factorizable state yields J-r{A, B) = 0. 

Second, we define geometric and arithmetic averages in 
order to obtain more information about the state and its 
type of entanglement. We define the arithmetic average 

M = N{V{A,B)){C^)-'Y.V{A,B) , (12) 



and the geometric average 

g=N{V{A,B)){j{v{A,B)y " , (13) 

where V{A,B) are the probe quantities: Qc{pab) and 
J^r{pAB) in our case, and the sum and product are over 
all possible pairs of qubits. We note that V{A,B) = 
V{B, A) and that the normalization factor, M{V{A, B)), 
is introduced so the measure yields one for generalized 
GHZ states. We have TV = (1 + (d - 1)(1 - S2.N)) if we 
are averaging J-'r or M — 1 if we are using the quasi- 
concurrence. 

The generalization to the case of mixed states can be 
achieved through the expressions Q 

M{p)=mmY^p,Mip,) , (14) 

and 

g{p) = mm^piGip,) , (15) 

for the arithmetical and geometrical averages respec- 
tively. The minimum is intended over all possible decom- 
positions. Note that in the two-qubit case, A4{p) = G{p) 
and they reduce to the concurrence defined by Woofers 
p^ . It is important to note that the probe quantities 
are directly applied without minimizations when they are 
acting as probes on a multipartite quantum state; when 
analyzing bipartite states, however, they act as full mea- 
sures so they operate through minimizing mechanisms 
such as the ones given above. 

Both quantities give different information about the 
type of entanglement of the quantum state, namely if 
it has global entanglement (among all parties) or only 
among some of the parties (bipartitions, tripartitions, 
etc.), and, in conjunction, they give us the possibility of 
fully characterizing the entanglement of a pure quantum 
state. 

The arithmetic average varies between zero and one: it 
equals zero if and only if all QciPABY^ (or J-r{A, B)'s) 
are zero, that is, if the state is completely factorizable, 
and equals one if and only if all Qc(pab)'s are 1 (or 
J^r{A, Bys are 1/2), i.e. if we have an 7V-qubit max- 
imally entangled state. It can be shown [20[ that the 
measure defined in this way is both additive and strongly 
super additive, two properties that are desired but usu- 
ally not satisfied by most entanglement measures [2]| . 
Thus, we provide a measure of entanglement that accu- 
rately quantifies the amount of entanglement of a pure 
quantum state. This measure is on its own, however, 
unable to determine whether a state is genuine globally 
entangled or not. For example, a bi-factorizable iV-qubit 
state could yield the same value as a globally entangled 
state. 

The geometric average also varies between zero and 
one. This equals one if and only if all Qc{pab)'s 
{Tr{A, B)^s) are 1 (1/2), that is, if we have a maximally 
entangled iV-qubit state. However, and in contrast with 
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the arithmetic measure, it equals zero if at least one of 
the 'P{A, B) vanishes, i.e., if the state is at least bifactor- 
izable, thus quantifying global entanglement. Higher fac- 
torizabilities, e.g. trifactorizable or 7V-factorizable states 
would have even more vanishing 7^ (A, i3)'s. 

These two measures lead us to classify the set of pure 
multipartite states between the ones with non vanishing 
Q, i.e. genuine globally entangled statesjl^, and the rest, 
which may or not be entangled at all, information that 
can be obtained through the arithmetical average M. . We 
recall that according to our definition a state is genuine 
globally entangled if after measuring one qubit we gain 
some information about all of the other qubits in the reg- 
ister, that is iiV{A, B) ^ for all {A, B). Both measures 
provide a faithful way of discriminating between pure 
multipartite quantum states. We can further explore the 
structure of the globally entangled states; there are two 
distinct types: the ones with the same 7^(A, B) for all 
pairs of qubits (homogeneously entangled states) and the 
states that posses different values of non-factorizability 
for all their pairs (heterogeneously entangled states). In 
the context of this classification, our calculations lead us 
to define a GHZ state as the homogeneously entangled 
state with the highest possible average of 7^(A, B). 

We introduce this definition because, if we consider 
the case of the quantity J^r j28| , given an arbitrary state 
we may find higher values of correlations between a pair 
of qubits, e.g., i) \EPR) (g) \EPKl states, and ii) max- 
imally entangled mixed states [22], but this does not 
mean that their degree of entanglement is higher than 
that of the \GHZ) state. For the case i) Tr{A,B) = 
1 > 1/2 for {A,B) = (1,2), (3, 4), and J^r{A,B) = 
for {A,B) = (1,3), (1,4), (2, 3), (2, 4), yielding M = 
2/3 < 1 = Mghz- Note that = < 1 = Gghz- 



For the case ii) consider the pure state \MEMS) = 
10101) + Vx/2(|0000) + 10011) + |1100) + |1111)), 
built as a purification of the maximally mixed entangled 
state proposed by Munro et al. (2^ : a direct calculation 
yields .Fr(l,2) = .Fr(3,4) > 1/2 = J^r{l,2)ipGHz) for 
X — > 1. Despite this result, it is not difhcult to check 
that J^r(l,3) = J^r(l,4) = J^r(2,3) = J^r(2,4) < 1/4 
which gives M < Mghz, and also Q < Qghz- In this 
way we see that, although these states exhibit a higher 
degree of non-factorizability among some of their pairs, 
on average they possess less entanglement than the GHZ 
state, which possesses the maximum possible entangle- 
ment on average. It is interesting that this result implies 
the possibility that weaker (less correlated) links may be 
taking place in spin chains or related systems [25j . 

We presented a geometrical formalism which is suited 
to describing the Meyer- Wallach [6| measure in its Bren- 
nen's (23j version, thus generalizing previous work by 
Mosseri 



17|. We also in- 



Bernevig [15[, Mosseri [16|, and Levay 
troduced, through the concept of probe quantities char- 
acterizing the non-factorizability of a bipartite density 
matrix, the construction of arithmetic and geometric en- 
tanglement measures for quantifying multipartite pure 
states, which were able to distinguish between globally 
and partially entangled states and that accurately quan- 
tified the degree of entanglement of an A^-qubit system. 

Acknowledgments. We thank T. Brandes, C. 
Emary, and H. Ocampo for useful discussions. JHR 
thanks S. Reina-Steers for useful distractions. GAPS 
thanks G. H. Paz, I. Silva, D. F. Gutierrez, and G. R. 
Paz for continual support. We gratefully acknowledge 
financial support from COLCIENCIAS under Research 
Grants No. 1106-14-17903 and No. 1106-05-13828. 



C. H. Bennett and D. P. DiVincenzo, Nature 404, 247 [17] 

(2000) . [18] 
See, e.g. M. A. Nielsen and I. L. Chuang, Quantum Com- [19] 
putation and Quantum Information, CUP (2000). 

C. H. Bennett et al, Phys. Rev. Lett. 70, 1895 (1993). [20] 
A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

M. B. Plenio an d S. Virmani, E-print [21] 
|arXiv:quant-ph/0504163 | 

D. A. Meyer and N. Wallach, J. Math. Phys. 43, 4273 [22] 
(2002). 

J. G. Luque and J. Y. Thibon, J. Phys. A 39, 371 (2006). [23] 
C. Emary, J. Phys. A 37, 8293 (2004). [24] 
V. Vedral et al, Phys. Rev. Lett. 78, 2275 (1997). [25] 

A. Peres, Phys. Rev. Lett. 77, 1413 (1996). [26] 
M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 

Lett. A 223, 1 (1996). [27] 

P. Levay, J. Phys. A 37, 1821 (2004); 38, 9075 (2005). 

G. Jaeger et al. Phys. Rev. A 68, 022318 (2003). [28] 

V. Coffman, J. Kundu, and W. K. Wootters, Phys. Rev. 

A 61, 052306 (2000). 

B. A. Bernevig and H. Chen, J. Phys. A 36, 8325 (2003). 
R. Mosseri and R. Dandoloff, J. Phys. A 34, 10243 

(2001) . 



5194 



P. Levay, J. Phys. A: Math. Gen. 37, 1821 (2004). 
A. J. Scott, Phys. Rev. A 69, 052330 (2004). 
N. J. Cerf and C. Adami, Phys. Rev. Lett. 79, 
(1997). 

G. A. Paz-Silva and J. H. Reina, E-print 
arXiv:quant-ph/0610017 

M. Christandl, D.Phil, thesis. University of Cambridge, 



E-print |arXiv:qua nt-ph/06 04183 

W. J. Munro, D. F. V. James, A. G. White and 
P. G. Kwiat, Phys. Rev. A 64, 030302 (2001). 
G. K. Brennen, Quant. Inf. and Comp. 3, 616 (2003). 
P. J. Love et al., E-print arXiv:quant-ph/0602143 
S. Bose, E-print arXiv:cond-mat/0610024J 
However, an invariant common to both can be con- 
structed using the norm of the vector. 
For higher A'^ we would have different partitions, e.g. 6 — 
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definition is even more compact; it would then read: a 
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